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Additive-State-Decomposition-Based Output Feedback Tracking Control 
for Systems with Measurable Nonlinearities and Unknown Disturbances 

Quan Quan, Kai-Yuan Cai, Hai Lin 



Abstract 

In this paper, a new control scheme, called additive-state-decomposition-based tracking control, 
is proposed to solve the output feedback tracking problem for a class of systems with measurable 
nonlinearities and unknown disturbances. By the additive state decomposition, the output feedback 
tracking task for the considered nonlinear system is decomposed into three independent subtasks: a 
pure tracking subtask for a linear time invariant (LTI) system, a pure rejection subtask for another 
LTI system and a stabilization subtask for a nonlinear system. By benefiting from the decomposition, 
the proposed additive-state-decomposition-based tracking control scheme i) can give a potential way 
to avoid conflict among tracking performance, rejection performance and robustness, and ii) can 
mix both designs in time domain and frequency domain for one controller design. To demonstrate 
the effectiveness, the output feedback tracking problem for a single-link robot arm subject to a 
sinusoidal or a general disturbance is solved respectively, where the transfer function method for 
tracking and rejection and the feedback linearization method for stabilization are applied together 
to the design. 



Index Terms 

Additive state decomposition, output feedback, measurable nonlinearities, tracking, rejection. 



I. Introduction 

In this paper, the output feedback tracking problem for a class of systems with measurable 
nonlinearities and unknown disturbances is considered. This problem has attracted great 
research interest in recent years [QQl-Gl- As far as nonlinear systems are concerned, several 
results are available under the minimum phase assumption. In [1], global disturbance rejection 
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with stabilization for nonlinear systems in output feedback form was solved in spite of the 
disturbance generated by a finite dimensional exosystem. The similar problem but subject to 
unknown parameters on both input matrix and system matrix was considered in B3. In 0, 
another control algorithm in simplicity of implementation was proposed for nonlinear plants 
with parametric and functional uncertainty in the presence of biased harmonic disturbance. 
The problem about nonminimum phase nonlinear systems was further considered. In [3]|- 
JH, adaptive estimation of unknown disturbances in a class of nonminimum phase nonlinear 
systems, and the stabilization and disturbance rejection based on the estimated disturbances 
for single-input single-output (SISO) systems were solved. The result was further extended 
to a class of nonminimum phase nonlinear Multiple-Input Multiple-Output (MIMO) systems 
in 0. In |5), a solution to this problem was provided for nonminimum phase systems with 
uncertainties in both parameters and order of an exosystem. 

The basic idea of the current work is to decompose the output feedback tracking task 
into simpler subtasks. Then one can design a controller for each subtask respectively, which 
are finally integrated together to achieve the original control task. The motivation of this 
paper can be described as follows. First, it is to avoid conflict among tracking performance, 
rejection performance and robustness. It is well known that there is an intrinsic conflict 
between performance (trajectory tracking and disturbance rejection) and robustness in the 
standard feedback framework [8J,[9|. By the control scheme mentioned in Q]|-[|7]|, as the 
dimension of the exosystem is increasing, the closed-loop system will be incorporated into 
a copy of marginally stable exosystem according to internal model principle IfTOl to achieve 
high performance (trajectory tracking and disturbance rejection). The price to be paid is a 
reduced robustness against uncertainties, especially for nonminimum phase systems. More- 
over, conflict between tracking performance and rejection performance exists as well when 
reference and disturbance behave differently [fTTTl . Secondly, it is to relax the restriction on 
the disturbance. For the control scheme mentioned in flU-®, a "regulator equation" often 
needs to be solved first for a coordinate transformation, which yields an error system with 
disturbances appearing at the input. However, these control schemes are only applicable 
to finite-dimensional autonomous exosy stems. While, the control scheme mentioned in Q 
requires the system being minimum phase to shift disturbances to the input channel. 
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For such a purpose, a new control scheme based on the additive state decomposition^, called 
additive-state-decomposition-based tracking control, is proposed which is applicable to both 
minimum phase and nonminimum phase systems. The proposed additive state decomposition 
is a new decomposition manner different from the lower-order subsystem decomposition 
methods existing in the literature, see e.g. [fT3ll . lfT4l . Concretely, taking the system x(t) = 
f (t,x) ,x G M. n for example, it is decomposed into two subsystems: x\ (t) = fi (t,xi,x 2 ) 
and ±2 (t) = f2 (t, xi, X2), where X\ G W lx and x 2 G IR™ 2 , respectively. The lower-order 
subsystem decomposition satisfies 

n = rii + 72,2 and x = x\ © x 2 - 

By contrast, the proposed additive state decomposition satisfies 

n = rii = n 2 an d x = X\ + x 2 - 

In our opinion, lower-order subsystem decomposition aims to reduce the complexity of 
the system itself, while the additive state decomposition emphasizes the reduction of the 
complexity of tasks for the system. 

By following the philosophy above, in the additive-state-decomposition-based tracking con- 
trol scheme, the output feedback tracking is 'additively' decomposed into three independent 
subtasks, namely the tracking subtask, the rejection subtask and the stabilization subtask. 
Three subcontrollers for the three subtasks are designed separately then. Since the resulting 
controller possesses three degrees of freedom, the proposed scheme in fact gives a potential 
way to avoid conflict among tracking performance, rejection performance and robustness. 
Moreover, by the additive-state-decomposition-based tracking control scheme, it will be seen 
that both the tracking subtask and rejection subtask only need to be achieved on a linear time 
invariant (LTI) system. Consequently, the tracking controller and disturbance compensator 
can be designed in both time domain and frequency domain. In this framework, the existing 
output regulation methods as in Q]|-[[6]| can be incorporated. Also, it can take advantage of 
some standard design methods in frequency domain to handle general disturbances. More 
importantly, nonminimum phase systems can be handled in the same framework. 

This paper is organized as follows. In Section 2, the problem formulation is given and 
the additive state decomposition is recalled briefly first. In Section 3, the considered system 

'in this paper we have replaced the term "additive decomposition" in 1121 with the more descriptive term "additive state 
decomposition". 
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is transformed to a disturbance-free system in sense of input-output equivalence. Sequently, 
in Section 4, the transformed system is 'additively' decomposed into three subsystems. In 
Section 5, controller design is given. Section 6 concludes this paper. 

II. Problem Formulation and Additive Decomposition 
A. Problem Formulation 

Consider a class of SISO nonlinear systems similar to (D-0|,|l6l-|ED: 

x = A x + bu + 4>o (y) + d,x (0) = xq 

V = c T x (1) 

where A E R nxn is a constant matrix, &GK" and ceR" are constant vectors, <p : R — > R n 
is a nonlinear function vector, x (t) & R n is the state vector, y (t) G R is the output, u (t) & R 
is the control, and d (t) e R n is a bounded disturbance. It is assumed that only y is available 
from measurement. The desired trajectory r(t) el is known and smooth enough, t > 0. In 
the following, for convenience, we will omit the variable t except when necessary. 
Remark 1. Under certain conditions, the system in the form 

x — f (x) + g (x) u + d 
y = h(x) 

can be transformed to (OQ). The sufficient and necessary condition to ensure the existence of 
transformation can be found in [fT5ll . 

Remark 2. The considered SISO nonlinear system (QQ) is allowed to be a nonminimum 
phase system, i.e., the transfer function of the linear part (i.e., regardless of nonlinear dy- 
namics 4> (y) and disturbance d) 

is nonminimum phase here, where iV (s) has zeros on the right s-plane. It is noticed that the 
property of nonminimum phase cannot be changed by output feedback. 

For system (0Q), the following assumption is made. 

Assumption 1. The pair (A ,c) is observable. 

Under Assumption 1, the objective here is to design a tracking controller u such that y — > r 
as t — > oo or with good tracking accuracy, i.e, y — r is ultimately bounded by a small value. 
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B. Additive State Decomposition 

In order to make the paper self-contained, additive state decomposition lfT2l is recalled 
here briefly. Consider the following 'original' system: 

f(t,x,x) = 0,x(0) = x (2) 

where x G M n . We first bring in a 'primary' system having the same dimension as ©, 
according to: 

fp (tj &pi %p) 0, Xp (0) %pfi (3) 

where x p G M n . From the original system © and the primary system © we derive the 
following 'secondary' system: 

/ (t, x, x) - f p (t, x p , x p ) =0,x (0) = x (4) 

where x p G IR n is given by the primary system ©. Define a new variable x s G M n as follows: 

Then the secondary system © can be further written as follows: 

/ (t, x s + Xp, x s + x p ) - f p (t, x p , x p ) = 0, x s (0) =x - Xp fi . (6) 

From the definition ©, we have 

x (t) = x p (t) + x s (t),t> 0. (7) 

Remark 3. By the additive state decomposition, the system © is decomposed into two 
subsystems with the same dimension as the original system. In this sense our decomposition 
is "additive". In addition, this decomposition is with respect to state. So, we call it "additive 
state decomposition". 

As a special case of ©, a class of differential dynamic systems is considered as follows: 

x = f(t,x),x(0) = x , 

y = h(t,x) (8) 

where x G W l and y G R m . Two systems, denoted by the primary system and (derived) 
secondary system respectively, are defined as follows: 

Xp fp (t, Xp) , Xp (0) %p t o 

V v = hp (t, Xp) (9) 
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and 

X s f (t, Xp -\- X s ) fp (t : Xp^j , X s (0) Xq 2^,0) 

y s = h (t, x p + x s ) - hp (t, x p ) (10) 

where x s = x — x p and y s = y — y p . The secondary system (flOl) is determined by the original 
system ([8]) and the primary system ©. From the definition, we have 

x (t) = x p (t) + x s (t) ,y(t)= y p (t) +y s (t),t> 0. (11) 

III. Model Transformation 

Firstly, we need to estimate the state from the output. The main difficulty is how to handle 
the disturbances in the state equation. In Q]|-|I2L an extended state observer including states 
of the considered nonlinear system and the exosystem is designed, where the disturbance 
is assumed to be generated by a finite-dimensional autonomous exosystem. According to 
this, the model of the disturbance has in fact determined the performance of the observation 
partly. However, in practice, a general disturbance is difficult to model as a finite-dimensional 
autonomous one, or with uncertainties. To tackle this difficulty, we first transform the original 
system (Q~|) to a disturbance-free system, which is proved to be input-output equivalent with 
the aid of the additive state decomposition as stated in the following theorem. 

Theorem 1. Under Assumption 1, there exists a vector p E W 1 such that A = A + pc T is 
stable, and the system (QQ) is input-output equivalent to the following system: 

inew = Ax new + bu + (j) (y) , x new (0) = 

y c x new -\- d new (12) 

where (y) = </> (y) — py and 

Xnew = x - (e At x + J e A{t ~ s) d (s) ds\ , d new = c T (e At x + J e A(t ~ s) d (s) ds \ . (13) 

Proof. Since the pair (A , c) is observable, there always exists a vector pel" such that 
A = A Q + pc T is stable, whose the eigenvalues can be assigned freely. The system (OQ) can 
be rewritten as follows: 

x = Ax + bu + 4> (y) + d,x (0) = x 

y = c T x (14) 
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where 4> (v) — 4>o iv) ~ PV- I n the following, additive state decomposition is utilized to 
decompose the system (fl4l) . Consider the system (PT41) as the original system and choose the 
primary system as follows: 

Xp = Axp ~\~ d^ Xp (0) = Xq 

y p = c T x p . (15) 



Then the secondary system is determined by the original system (|T4l) and the primary system 
(TBI) with the rule dJO]) that 

x s = Ax s +bu + (j){y) , x s (0) = 

y s = c T x s . (16) 

According to (fTTI) . we have x = x p + x s and y = y p + y s - Therefore, the following system 
is an input-output equivalent system of (fl4l) : 

x s = Ax s + bu + <p(y) , x s (0) = 

y = c T x s + c T x p (17) 

where x p is generated by (fl"5T) . By (fl"5l) . it holds that x p = e At x + L e j4 ^~ ,s ^ (s) <is. Let 
= %new and <i ne ^ = c T x p . Substituting into (TTTI) yields (fT2l) . □ 

For the disturbance-free transformed system (fT2l) . we design an observer to estimate x ne „, 
and d new , which is stated in Theorem 2. 

Theorem 2. Under Assumption 1, an observer is designed to estimate state x new and d new 
in (fT2l) as follows 

%new ^^'Jiety "I - (?/) ) %new (0) 

dnew 2/ C %new (18) 



Then x new — •Enew and d new — d new . 

ZVoo/ Subtracting (IT8T) from (fT2l) results in x ne „, = Ax new , x new (0) = 0, where x nCT() = 
x neu) — x new . Then x neu , = 0. This implies that x new = x new . Consequently, by the relation 
y = c T x new + d new in dT2j), we have d new = d new . □ 

Remark 4. By (fT3T) . if the new state x new is bounded, then the original state x is bounded 
as well since the matrix A is stable and the disturbance d is bounded. This explains why the 
matrix A is chosen to be stable. To eliminate the transient effect of initial values in d new , 
namely c T e At x , we often assign the eigenvalues for A to have large negative real part. By 
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using the new state x new , the controller can be designed based on the transformed system 
(fl"2l) directly as shown in the following sections. 

Remark 5. It is interesting to note that the new state x new and disturbance d new in the 
transformed system (fl"2l) can be observed directly rather than asymptotically or exponentially. 
This will facilitate the analysis and design later. In practice, the output y will be more or less 
subject to noise. In this case, the stable matrix A will result in a small x new in the presence 
of small noise, i.e. x new close to x new . 

Example 1. A single-link robot arm with a revolute elastic joint rotating in a vertical plane 
is served as an application in this paper IfToll : 



X\ = x 2 

Ft k Mgl . 

x 2 = — -x 2 (xi — x 3 ) — sin x\ + d\ 

Ji Ji Ji 

x 3 = 

F m k 1 
£4 = — + — (xi - x 3 ) + —t + d 2 

Jm Jm Jm 

y = x l (19) 

where xi,x 2 ,x 3 ,x 4: are the link displacement (rad), link velocity (rad/s), rotor displacement 
(rad) and rotor velocity (rad/s), respectively; d 1 and d 2 are unknown disturbance. The initial 

1 T 



value is assumed to be x (0) 



0.05 0.05 



. Let link inertia J\ = 2kg -m , the 



motor rotor inertia J m = 0.5kg-m 2 , the elastic constant k = 0.05kg-m 2 /s, the link mass 
M = 0.5kg, the gravity constant g = 9.8m/s 2 , the center of mass I = 0.5m and viscous 
friction coefficients Fi = F m = 0.2kg-m 2 /s. The control r is the torque delivered by the 
motor. The control problem here is: assuming only y is measured, r is to be designed so that 
y tracks a smooth enough reference r asymptotically or with good tracking accuracy. The 
controller r in (fT9l) is designed as follows r = J m u, where u will be specified later. Then 
the system (fT9l can be written in form of © with 








1 















1 







A = 


k 
Ji 




Fi 
Ji 




k 

Jl 






1 


,b = 






,c = 






, 0o (y) = 


-M^siny 





k 

Jm 





k 


F m 

Jm 




1 




_ _ 








(20) 

It is easy to verify that the pair (A , c) is observable. So for this application Assumption 1 

T 



holds. It is found that A is unstable. Choosing p 



-2.10 -1.295 -9.36 3.044 
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the system (1191) is formulated into (1121) with 

-2.1 10 

-1.32 -0.1 0.025 

-9.36 1 

3.144 -0.1 -0.4 



.4 



(y) 



2.1y 

1.2952/ — 1-225 sin y 
9.36y 
-3.044^/ 



where the eigenvalues of A are assigned as —0.5, —0.6, —0.7, —0.8. 



IV. Additive State Decomposition of Transformed System 



In this section, the transformed system (1121) is 'additively' decomposed into three in- 
dependent subsystems in charge of corresponding subtasks, namely the tracking subtask, 
the rejection subtask and the stabilization subtask, as shown in Fig.l. There exist many 
tools to analyze LTI systems, such as Laplace transformation and transfer function, and 
state-space techniques. Based on the above consideration, the transformed system (PT2l) is 
expected to be decomposed into two subsystems by the additive state decomposition: an 
LTI system including all external signals as the primary system, together with the secondary 
system free of external signals. Therefore, the original tracking task for the system (fT2l) is 
correspondingly decomposed into two subtasks by the additive state decomposition: a tracking 
(including rejection) subtask for an LTI 'primary' system and a stabilization subtask for the 
left 'secondary' system. Since the tracking (including rejection) subtask is only assigned to 
the LTI system, it therefore is a lot easier than that for a nonlinear one. Furthermore, the 
tracking (including rejection) subtask is decomposed into a pure tracking subtask and a pure 
rejection subtask. 

Consider the transformed system (PT21) as the original system. According to the principle 
above, we choose the primary system as follows: 



x p = Ax p + bu p + (j) (r) 

Up c Xp -\- d newi Xp (0) 0. 



(21) 



Then the secondary system is determined by the original system (fl2l) and the primary system 
(I2TI) with the rule (flOl) . and we can obtain that 



Ax s + bu s + (j) (r + e p + c T x s ) - (r) 
c T x s ,x s (0) = 



(22) 
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Tracking (Rejection) Task on 
Nonlinear System (12): y —> r 



£ 



Tracking (Rejection) Subtask 
on LTI System (21): y -> r 



£ 



y p y pp y ps 



Tracking Subtask on LTI 
System (25): y -> r 



y = y P + y s 



1 



Stabilization Subtask on 
Nonlinear System (22): y s -> 



Rejection Task on LTI 
System (26): y ps -> 



Fig. 1. Additive state decomposition flow 



where u s = u — u p and e p = y p — r. According to (fTTI) . we have 

Xnew = Xp + x s and y = y p + y s . (23) 

The strategy here is to assign the tracking (including rejection) subtask to the primary 
system (|2TT) and the stabilization subtask to the secondary system (l22l) . It is clear from (f2Th - 
d23l that if the controller u p drives y v — )■ r and the controller u s drives y s — > as t — > oo, 
then y — > r as t — > oo. The benefit brought by the additive state decomposition is that the 
controller u s will not affect the tracking and rejection performance since the primary system 
(I2TI) is independent of the secondary system (1221) . On the other hand, if the secondary system 
(1221) is input-to-state stable with respect to the signal e p , then the dynamics of controller u p 
will not change the input-to- state stability property. Therefore, conflict between performance 
(trajectory tracking and disturbance rejection) and robustness is avoided. Since the states 
x p and x s are unknown except for sum of them, namely x new , an observer is proposed to 
estimate x p and x s . 

Theorem 3. Under Assumption 1, suppose that an observer is designed to estimate state 
x p and x s in (l2"TI)-(l22l) as follows: 

X p %new x s (24a) 
x s = Ax s + bu s + (f) (y) - cf) (r) , x s (0) = 0. (24b) 

Then x p = x p and x s = x s . 
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Proof. Subtracting (I24bl) from (|22l) results in x s = Ax s , x s (0) = [1 where x s = x s — x s . 
Then x s = 0. This implies that x s = x s . Consequently, by (1231) . we have x p = x new — x s = x p . 

□ 

To avoid conflict between tracking performance and rejection performance, the tracking 
subtask for the primary system (12TI) is further decomposed into two subtasks by the additive 
state decomposition: a pure tracking subtask and a pure rejection subtask. At this time, 
consider the primary system (TfTI) as the original system and choose the primary system of 
(f21~b as follows: 

Xpp A.Xpp ~\- bzipp ~\- (j) (/*) ; 

Vpp c ■Eppi Xpp (0) 0. (25) 

Then the secondary system of (12TI) is determined by the system (12TT) and (1251) with the rule 
([TO]) that 

j\xp S -\- \yiip S 

Vps c -Eps d newi Xp S (0) 0, (26) 

where w ps = u p — u pp . According to (fTTT) . we have 

Xp -|- Xp S and t/p ?/pp -|- Up S . (27) 

It is clear from (l25l)-(l271) that if the controller u pp drives y pp — > r and the controller u ps drives 
y P s — > as t — > oo, then y p — > r as t — > oo. It is noticed that the controller w P p and u ps 
above are independent each other. So, conflict between tracking performance and rejection 
performance is avoided. 

V. Controller Design 

So far, we have transformed the original system to a disturbance-free system whose state 
can be estimated directly. And then, decompose the transformed system into three independent 
subsystems in charge of corresponding subtasks. In this section, we are going to investigate 
the controller design with respect to the three decomposed subtasks respectively. 

2 Since the initial values x new (0) , x a (0) , x s (0) are all assigned by the designer, they are all determinate. 
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A. Problem for Tracking Subtask 

Problem 1. For rf25l) . design a control input 

Z Oi (z , Xppj Tj 

u pp = u r (z r ,x pp ,r) (28) 

such that y pp — > r as t — > oo. 

Remark 6 (on Problem 1). Problem 1 can be considered as a stable inversion problem 
lfT7Tl . lfT8ll . which has been solved for the system in the form of (|25l) no matter whether it 
is minimum phase or nonminimum phase. If the system (1251) is only minimum phase, then 
the stable inversion problem is made a lot easier by using the transfer function method. 
Problem 1 can also be considered as an output regulation problem Q]|-[[6]|, if the reference 
r is generated by an autonomous system. Since the reference r is given, the computation 
above can be done offline. In addition, since the control design is only for a classical LTI 
system, the computation is easier compared with that for a nonlinear one. 

Example 2 (Example 1 Continued). Since system (TT9l) is minimum phase, by the transfer 
function method, the reference control input u pp = u r can be designed as follows: 



u r (s) 



r (s) + c T F 1 (s) ( pr (s) + b\ — j— (sin r) (s 

Ji 



(29) 

T 



where u r (s) is the transfer function of signal u r , F (s) = sI^ — A and foi=^0 1 oj- 
Then y pp — > r as t — > oo. Problem 1 is solved. 

B. Problem for Rejection Subtask 

Problem 2. For ft26\l . there exists a control input 

z d = a d (z d r rl \ 

U ps U i^Z , X pS) dnew) (30) 

such that y ps — > B (5) □ as t — >■ oo. In particular, if 5 = 0, then y ps — > as t — > oo. 

Remark 7 (on Problem 2). Since the system (T26l) is a classical LTI system, some standard 
designs in frequency domain, such as the transfer function method, can be used to handle a 
general disturbance IfTTTl . In this case, the disturbance cannot be often rejected asymptotically. 

3 B(S) = £ ^ 8}} where S — S (d) is a function of the disturbance d; the notation x (t) — > 13(6) means 

min \x (t) — y \ — s> 0. 
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So, the Problem 2 needs to consider the result y ps — > $ (5) besides j/ ps — >■ 0. If d new is 
generated by an autonomous system, then Problem 2 can be considered as an output regulation 
problem flTJ -[JVj . In this case, the disturbance can be rejected asymptotically. The technique 
in ttTJ-0 of course can be still applied to the problem even if both parameters and order of 
an exosystem are uncertain. 

Example 3 (Example 1 Continued). To demonstrate the effectiveness of the proposed con- 
trol, the disturbances in ([191) are assumed to be in two cases. Case 1 (sinusoidal): the unknown 
disturbances d\ and d 2 are sinusoidal, such as d\ = 0.05 sin t and d 2 = 0.05 sin (t + 1) . 
Case 2 (general): the unknown disturbances d\ and d 2 are driven by normal (Gaussian) 
distributed random signals (1,(2 ~ jV(0, 1). The transfer functions are assumed to be 



d x is) 



0.08 



rCi (s) and d 2 (s) 



6 



(s+l)(s+2)(s+3) 



( 2 (s). The resulting disturbances 



(s+0.2)(s+0.5)(s+0.8) 

are shown in Fig. 2. 

Case 1 (sinusoidal). Since the unknown disturbances d\ and d 2 are sinusoidal with fre- 
quency 1 rad/s, by (PT31) we can conclude that d new = d + e, where d is sinusoidal with 
frequency 1 rad/s as well and e — > as t — > 00. This implies that the transfer function d can 
be written as d(s) = -^iw (0) . For (|26l) . design controller follows 



Sz d + c d ?/p S 



yd 1 i,^ "™. 



ps 



where S = 
Combining 



1 

-1 
with 







96 122 



T 



and k 2 



(3D 



34 70 -0.67 -1.9 



3TI) results in 



+ CdC x ps + c d d. 



X 



ps 



x ps = k[z d + (A + bkl) 
Ups c ■Eps d new . 

The transfer function from d new to y ps possesses negative real poles and at least two zeros 
±i. Since d(s) = gA^w (0) and e — > 0, we have y ps — > 0. Problem 2 is solved for Ca^e 1. 
Case 2 (general). The transfer function from u ps to y ps can be represented as follows: 



Vps \s) 



c T F~ l (s) bu ps (s) + d 



Since d new is a low-frequency disturbance and can be observed by (PT8l) . an easy way is to 
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design the disturbance compensator as follows: 

u d (s) = - — - rr Q(s)d new (s) (32) 
c 1 r ~ x [s) o 

where u d (s) is the transfer function of u d , Q (s) is a low-pass filter satisfying 

(1 - Q (ju)) d new (jlu) w 0, Vu; G R. 

Moreover, Q (s) is at least fourth order to make the compensator physically realizable (the 
order of denominator is greater than or equal to that of numerator). In this simulation, we 
choose 

Q 00 ' 



nsys'+i)' 

In this case, y ps in (1261) is ultimately bounded by a small value, namely y ps — >■ $ (<5), where 
5 can be adjusted by Q (s) . Problem 2 is solved for Case 2. 



C. Problem for Stabilization Subtask 

Problem 3. For ([22]), there exists a controller u s = u s (x s ,r, ■■■ , r^) such that the 
closed-loop system is input-to-state stable with respect to the input e p , namely 

\\x s {t)\\ </3 (\\x s (t )\\,t -t ) + 7 ( sup ||ep(s)|| ] ,t >t , (33) 

\t <s<t J 

where r^ N > denotes the Nth derivative of r, function (3 is a class K,C function and 7 is a 
class K function / f79l/ . 

Remark 8 {on Problem 3). If e p is nonvanishing, then Problem 3 is a classical input-to-state 
stability problem. Readers can refer to lfT9l . ll20ll for how to design a controller satisfying input- 
to-state stability or how to prove the designed controller satisfying input-to-state stability. In 
particular, if e p — > as t — > 00, then x s — > as t — > 00 by (1331) . In addition, if e p — > 
as t — > 00, then input-to-state stability can be relaxed as well. In fact, Problem 3 only 
considers how x s behaves as e p — > as t — > 00. The reference [12111 discussed under what 
conditions the solution of (1261) say xt F and the solution of (1261) with e p = say x* satisfy 
\\x e s p (t) — x* (t)\\ < 9e~ vt , where 9,r] > 0. In this case, only stability of (|26|) with e p = 
needs to be considered rather than input-to-state stability. 
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Example 4 (Example 1 Continued). The system (|22l) can be rewritten as 



X s ,l — Xs,2 

Fi 



X s ,2 



j-X s ,2 7" \ x s,l — Xs,3) 



— — [sin (r + x S)1 ) - sin (r)J + d ep 
Ji 



Xs,3 ^s,4 



X S A 



XsA + 



k 



(x St i - x S)3 ) + w s 



7 ' 7 

y = x s>1 (34) 

where <i ep = — [sin (r + e p + cc S) i) — sin (r + x a ,i)\ • By the feedback linearization method, 
design u s as follows 

J, 

(35) 



w s r, f, f) = /ii + — (v + // 2 ) 



where 



u = -7.5x Si i - 19a; Si2 - 17% - 77y 4 

= -V3 + ^-^s,l - -j~X s - 3 - -j- X s ,4 
'J-m -Jm Jm 

F i , M ^ / I -\ / , \ Mgl r 2 . .. . 2 . 

jU 2 = —-^4 H (??3 + r ) cos (x Sj i + r) — l{z s ,2 + ?") sm {x„ t i + rj + r cos (r) - r sin 

<A Ji Ji 

Fi k Mgl 

r)3 = — T x s,2 - -r {x s ,i ~ x s ,3) ;— [sm (x S! i + r) - sm (r)\ 

Ji Ji Ji 

Fi k Mgl 

r) 4 = — -7? 3 - — (x S: 2 - x sA ) — [(x s>2 + r) cos (x s>1 + r) - r cos (r)J 

Ji Ji Ji 

r 1 T 

x s = x St \ x s ^ a; s ,3 x s,4 
Substituting (|35]> into ([34]) results in 

where 



x q x „ I dp 





X S ,1 







1 















X s ,2 










1 











,A' = 














V3 













1 













-7.5 


-19 


-17 


-7 








Since the matrix A' is stable and 
a class /C function 7' such that [[19 



— ~2fi W e pW' m ere exist a class K.C function j3' and 



<P'(\\x' s (t )\\,t-t )+i [ sup ||ep(a)|| ) ,t>t . 

y t <s<t 



Furthermore, by the definition of x' s , the Problem 3 is solved. 
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D. Controller Integration 

With the solutions of the three problems in hand, we can state 

Theorem 4. Under Assumption 1, suppose i) Problems 1-3 are solved; ii) the controller 
for system (QQ) (or (fl"2l) ) is designed as 
Observer: 

%new Ax new -\- blL ~\- (j) (?/) , X new (0) 

x s = Ax s + bu s + (y) - <p{r) , x s (0) = 0. 
x pp = Axpp + bu r + <fi (r) , x pp (0) = 



j 7 1 ^ ~ ~ a ~ ~ 

u neu; y ^ ^newi ^ p ^new lL ' s ; ^ps - L p ^ L, pp 



(37) 



Controller: 



z r = a r [z r , £pp, r) , 2 r (0) = 



z d 



u 



a 



Z i %ps ; d new i . 



z d (0) = 



^ i ^ti.s i d< 



Up + u s (£ s ,r, ■ ■ -r (Ar) ) . 



(38) 



Then the output of system (OQ) (or (fLU) ) satisfies that t/ — )■ r + £> (5 + ||c|| 7 (5)) as t — > 00. 
In particular, if 5 = 0, then the output in system (OQ) (or (fT2l) ) satisfies that y — > r as t — > 00. 
Proof. See Appendix. □ 

Remark 9. The controllers u pp and iip S tire designed based on the LTI systems, to which 
both design methods in frequency domain and time domain can be applied. The controller 
([37l)-d38l) has both following salient features: i) three degrees of freedom offered by three 
independent subcontrollers u r , u d and u s . This is similar to the idea of two-degree-of-freedom 
control ifTTI . If r and d new behave differently, then controller can be chosen both for good 
tracking of reference r and good rejection of disturbance d new . In addition, as the reference 
and/or disturbance change, only the corresponding subcontroller needs to be modified rather 
than the whole one. ii) The control signal u r and u d driven by reference and disturbance, 
considered as feedforward, will not effect on the stability of the closed-loop system, iii) The 
controller can deal with more general reference and disturbance signals more easily because 
subcontrollers u r and u d are designed based on a simple LTI system. 
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Example 5 (Examples 1-4 Continued). According to (1381) we design the controller as 
follows: 

z S z ~\~ Cd Xps ~t~ d neu ^ 



Case 1( sinusoidal): 



Case 2( general): 



Up — U -\- {k\ Z -\- 1^2 Xps) 

U = Up + u s 
Up = u r + u d 
U = Up + u s 



where x p , x s , x pp , d new are given by (l37l) with the parameters in (1201) . u r is given by (12~9l) . 
u s is given by (1331) and the disturbance compensator u d in Case 1 is designed according to 
(13TI) . while in Case 2 the disturbance compensator u d is designed according to (13~2l) . 

First, let us see the control performance of the primary system (l2~Tl) and the secondary 
system (1221) in Case i (sinusoidal disturbance). The evolutions of their outputs are shown 
in Fig. 3. As shown, it can be seen that y p = y pp + y ps — > r as t — > oo and y s — > as 
t — 7- oo. Therefore, by the additive state decomposition, y = y p + y s — > r as t — > oo. This is 
confirmed by the simulation shown in Fig.4. Secondly, let us see the control performance of 
the primary system (|2TI) and the secondary system (1221) in Case 2 (general disturbance). The 
evolutions of their outputs are shown in Fig. 5. As shown, it can be seen that y p tracks r with 
good performance and y s is ultimately bounded by a small value. Therefore, by the additive 
state decomposition, y = y p + y s tracks r with a good performance. This is confirmed by the 
simulation in Fig. 6. 

Remark 10. By the additive-state-decomposition-based tracking control scheme, it is seen 
from the simulation that the transfer function method is applied to the tracking controller 
design, which increases flexibility of the design. By benefiting from it, both a sinusoidal and 
a general disturbance can be handled in the same framework, where only the subcontroller 
u d needs to be modified rather than the whole one. 

VI. Conclusions 

In this paper, the output feedback tracking problem for a class of systems with measurable 
nonlinearities and unknown disturbances was considered. Our main contribution lies in the 
presentation of a new decomposition scheme, named additive state decomposition, which not 
only simplifies the controller design but also increases flexibility of the controller design. 

The proposed additive-state-decomposition-based tracking control scheme was adopted to 
solve the output feedback tracking problem. First, the considered system was transformed to 
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an input-output equivalent disturbance-free system. Then, by the additive state decomposition, 
the transformed system was decomposed into three subsystems in charge of three independent 
subtasks respectively: an LTI system in charge of a pure tracking subtask, another LTI system 
in charge of a pure rejection subtask and a nonlinear system in charge of a stabilization sub- 
task. Based on the decomposition, the subcontrollers corresponding to three subsystems were 
designed separately, which increases the flexibility of design. To demonstrate its effectiveness, 
the proposed additive-state-decomposition-based tracking control was applied to the output 
feedback tracking problem for a single-link robot arm with a revolute elastic joint rotating 
in a vertical plane. 

VII. Appendix: Proof of Theorem 4 
It is easy to follow the proof in Theorems 2-3 that the observer dTTT) will make 

&new = ^new, d n ew = d neW: Xp = Xp, 

/ >"■ , y - ry If If — — r f ( ^ ^ 

A s — ^pp — -^ppi ^ps — -"ps- \ Jy J 

The remainder proof is composed of two parts: i) for (I2TT) . the controller u p drives y p — > 
r + B (5) as t — > oo, and ii) based on the result of i), for (|22l) . the controller u s drives 
y s — >• B (||c|| 7 (5)) as t — > oo. Then the controller u = u p +u s drives y — > r+B (5 + ||c|| 7 (5)) 
as t —7- 00 in system (OQ) (or (TT21)). 

i) Suppose that Problems 1-2 are solved. By (|28l) and d39l) . the controller w pp is designed 
as follows: 

z r = a r (z r , x pp , r) , z r (0) = 

Wpp ^ ; ^ppj ?") 

which can drive y pp — > r as t — > 00 in (1251) . By (1301) and (|39l ), the controller -u ps is designed 
as follows: 

Z Oi (^Z , Xp S , dnew^j > ^ (0) 
tips ^ 3 &psi ^Tieiuj 

which will drive y ps — > i3 (5) as t — > oo in (1261) . Combining the two controllers u pp and u ps 
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above results in the controller for the primary system (|2~TT) : 

z r = a r [z r , x pp , r) , z r (0) = 

Z OL {^Z , Xp S1 dnewj ) % (0) 

"Up M (yZ , Xpp, 7") -|- "U ^2 , Xp S , f^new^ • (40) 

Therefore, by (1271) . the controller (1401) can drive y v — > r + B (5) as t — > oo. 

ii) Let us look at the secondary system (|22l) . Suppose that Problems 3 is solved. By (l39l) . 
the controller u s = u s (x s , r, • • • , r^) can drive the output y s such that 

lb. (*)|| < ||c|| ||x s (t)\\ 

< ||c||y3(||x s (to)||,*-*o) + ||c||7( sup \\e p (s)\\ j , t > 0. 

\to<s<* / 

Based on the result of i), we get e p — >■ i3 (5) as t — )■ oo. This implies that ||e p (t)\\ < 5 + e 
when t > t + Ti . Then 

< ||c|| /? (t + ,t — t - T\) + ||c|| ^ ( sup ||e p (s)|| ) ,t>t + T u 

< ||c|| /? (||x s . (t + TO || ,t - t - TO + ||c|| 7 (5 + e) , t > t + T x . 

Since ||c|| /3 (||x s (t + T\) || ,t — t — T x ) — ► as t — >■ oo and e can be chosen arbitrarily 
small, we can conclude y s — > B (||c|| 7 (5)) as t — )■ 00. Since y = c T x p + c T x s , we can 
conclude that, driven by the controller (1381) . the output of the system (OQ) (or (fT2l) ) satisfies 
that y — >■ r + B (5 + ||c|| 7 (5)) as t — )■ 00. In particular, if 6 = 0, then the output in system 
((O) satisfies that y — > r a.s t — > 00. □ 
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Fig. 2. Disturbance in Case 2 
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Fig. 3. Outputs of primary system and secondary system in Case 1 
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Tracking Performance 
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Fig. 4. Output (link displacement) of the single-link robot arm and control toque in Case 1 
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Fig. 5. Outputs of primary system and secondary system in Case 2 
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Tracking Performance 
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Fig. 6. Output (link displacement) of the single-link robot arm in Case 2 
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